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Abstract 

It is suggested that the propagation of the action potential is accompanied by an axoplasmic 
pressure pulse propagating in the axoplasm along the axon length. The pressure pulse stretch- 
modulates voltage-gated Na + (Nav) channels embedded in the axon membrane, causing their 
accelerated activation and inactivation and increasing peak channel conductance. As a result, the 
action potential propagates due to mechano-electrical activation of Nav channels by straggling 
ionic currents and the axoplasmic pressure pulse. The velocity of such propagation is higher than 
in the classical purely electrical Nav activation mechanism, and it may be close to the velocity of 
propagation of pressure pulses in the axoplasm. Extracellular Ca 2+ ions influxing during the voltage 
spike, or Ca 2+ ions released from intracellular stores, may trigger a mechanism that generates 
and augments the pressure pulse, thus opposing its viscous decay. The model can potentially 
explain a number of phenomena that are not contained within the purely electrical Hodgkin- 
Huxley-type framework: the Meyer- Overton rule for the effectiveness of anesthetics, as well as 
various mechanical, optical and thermodynamic phenomena accompanying the action potential. It 
is shown that the velocity of propagation of axoplasmic pressure pulses is close to the measured 
velocity of the nerve impulse, both in absolute magnitude and in dependence on axon diameter, 
degree of myelination and temperature. 
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I. INTRODUCTION 



The well-known Hodgkin-Huxley (H-H) model of propagation of the nerve impulse in 
neuronal axons [l| has recently been the subject of constructive criticism {2, 3|. As a purely 
electrical model, the H-H model does not contain a number of non-electrical phenomena 
which need to be included in a more complete model of axon function. These are changes 
in nerve dimensions and in the normal force exerted by the nerve 4|-[8|, reversible changes 
in temperature and heat 9[_ll3 1 and changes in fluorescence intensity and anisotropy of 
lipid membrane markers |l4j . Il5| . all observed synchronously with a propagating action 
potential. Interestingly, within the context of presently accepted models it has not been 



possible to explain the famous Meyer- Overton rule |!6l.ll7j for the effectiveness of anesthetics 
on nerve fibers. The rule states that the effectiveness of an anesthetic is linearly related to 
the solubility of that anesthetic in membranes 18) . The rule is valid over five orders of 
magnitude and holds independently of the chemical identity of the anesthetic; it is valid for 
noble gases such as argon and xenon as well as for alkanols. Although some ion channels are 
influenced by some anesthetics, there is no quantitative correlation with the Meyer-Overton 
rule J3], which makes it difficult to interpret the rule within the framework of the H-H type 
excitation kinetics. Although H-H type models have been very successful at reconstructing 

nn 

various aspects of action potential propagation H9J, new mechanisms may need to be 
considered to explain the Meyer-Overton rule. 

In this paper, we suggest an extension of the H-H model of the propagation of the 
action potential. We suggest that its propagation may be aided by stretch-modulation 
of membrane voltage-gated channels by an axoplasmic pressure pulse. This may lead to 
the action potential that propagates at a velocity close to the velocity of propagation of 
the axoplasmic pressure pulse. For both myelinated and unmyelinated fibers of various 
diameters, we estimate the pressure pulse velocity and indeed find that it agrees very well 
with experimental values for the velocity of the nerve impulse. The model offers explanation 
for the Meyer-Overton rule as well as for various mechanical, optical and thermodynamic 
phenomena that accompany the action potential. The speculation that density pulses may 



be related to nerve impulses has been discussed by various authors since 1912 20j-|23|. 



Most recently, it has been suggested 



3] that the action potential is a mechanical soliton 



(density pulse) propagating in the axon membrane. While this model appears to account 
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Propagating axoplasmic pressure pulse 
distends lipid bilayer, prior to activation 
of membrane Nav channels by straggling 
ionic currents 



Radial contraction of the lipid bilayer 
generates or augments propagating 
axoplasmic pressure pulse 
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Bilayer tension stretch-modulates Nav 
channels, enabling faster propagation of the 
depolarization wave in the stretched area 



Contracting actin-myosin filaments force 
a radial contraction in the lipid bilayer 
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Nav channels are activated by the combined 
action of straggling ionic currents and 
bilayer stretch 



Free intracellular Ca 2+ activates 
contraction in actin-myosin filaments 
anchored to the membrane interior 
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Ca 2+ is released from 
intracellular stores 



t 



A local H-H type voltage spike develops 



t 



The transmembrane voltage spike opens 
Cav channels, leading to Ca 2+ influx along 
the inward electrochemical gradient 



Nav channels are activated by depolarizing 
currents (e.g., at the axon hillock) 



FIG. 1: The suggested process of propagation of the action potential and the axoplasmic pressure 
pulse. 

for the above-mentioned phenomena, as well as for the Meyer-Overton rule, it does not seem 
to explain changes in the propagation velocity that are related to axon diameter. 



II. THE PRESSURE WAVE MODEL 



We suggest that the propagation of the nerve impulse may involve an axoplasmic pressure 
pulse propagating roughly synchronously with the wave of membrane depolarization. While 
several variations of the process can be envisioned, we suggest that it proceeds as outlined in 
Fig. 1. The axoplasmic pressure pulse is initiated by an action potential as discussed later 
in this section. The propagating axoplasmic pressure pulse distends the axon membrane 
prior to activation of membrane voltage-gated Na + (Nav) channels by the straggling ionic 
currents. This can occur if the axoplasmic pressure pulse propagates faster than the action 
potential associated with the H-H process for unstretched membrane. The stretch in the 
lipid bilayer modulates Nav channels, allowing a faster propagating H-H type excitation 
wave in the stretched membrane by enabling accelerated channel activation and inactivation 
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Radial bilayer 
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Axoplasmic pressure 
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FIG. 2: Amplification of the axoplasmic pressure pulse by a radial contraction of the axon mem- 
brane (drawn not to scale). 

and increasing peak channel conductance, as well as possibly lowering channel activation 
threshold. If this faster propagating H-H excitation wave is at least as fast as the traveling 
pressure wave, the two waves will co-propagate. In such a case, activation and development of 
local voltage spikes will be driven by the combined action of the membrane depolarization by 
the straggling ionic currents and the membrane stretch caused by the propagating pressure 
pulse. Further, we suggest the following mechanism of generation and amplification of the 
axoplasmic pressure pulse. (Below it is shown that under physiological conditions axoplasmic 
pressure pulses decay over roughly 1 mm distance due to viscosity, and therefore sustained 
propagation in long axons requires a mechanism of amplification.) A local H-H type voltage 
spike (which may be caused by an initial membrane depolarization such as at the axon 
hillock, or by the propagating action potential) activates membrane voltage-gated Ca 2+ 
(Cav) channels, leading to an influx of Ca 2+ ions along their inward electrochemical gradient; 
or, the voltage spike triggers Ca 2+ release from intracellular stores. The presence of free 
intracellular Ca 2+ induces contraction of actin-myosin filaments anchored to the membrane 
interior. The filament contraction forces a radial contraction in cylindrical segments of the 
membrane (Fig. 2), which generates or augments the propagating pressure pulse. 

Figure 3 illustrates suggested processes of mechanical action of the axoplasmic pressure 
pulse on Nav channels. The lipid bilayer is stretched laterally (Fig. 3(a)) as the propagating 
pressure pulse causes a AR(t) increase in the membrane cylinder radius R, corresponding to 
a Ar(i) increase in the lipid bilayer tension r. Tension r can be expressed through the bilayer 



lateral pressure profile p(z) 24j: r = f 2 h p(z)dz, where z is the depth within the bilayer, 



and the integral is over the bilayer thickness h. The change Ar(t) in the bilayer tension 
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to cell interior 



FIG. 3: Mechanical action of a traveling axoplasmic pressure pulse on Nav channels (see text). 

h 

results in a change Ap(z,t) in the lateral pressure profile such that Ar(t) = j\ Ap(z,t)dz. 
We suggest that these alterations in the bilayer lateral pressures modulate voltage-gating 
properties of Nav channels embedded in the bilayer. Alternatively, it is plausible that 
Nav channels are modulated by direct mechanical links (such as filaments) tethering the ion 
channel structures to the interior cytoskeleton that is perturbed by the propagating pressure 
pulse, as shown in Fig. 3(b). 

The model presented requires that Nav channels be mechanosensitive and allow a (per- 
haps substantially) faster propagating H-H excitation wave in a stretched membrane. Nav, 



Kv and Cav channels all exhibit reversible gating changes with stretch 



25 



261 ] . In experi- 



ments with pipette pressure applied to a membrane patch 26J], Nav channel activation and 
inactivation was reversibly accelerated by up to 1.85-fold (1.4-fold for typical stimuli used). 
Stretch also increased peak I^ a by about 1.5-fold near the activation threshold for moderate 
stretch stimuli. It is clear that such significant changes in Nav gating should lead to a faster, 
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and probably much faster, propagating H-H excitation. The above-mentioned experiments 
used stretch stimuli lasting on the order of tens of seconds. It is not excluded that for much 
shorter, ~1 ms, stretch stimuli, Nav channels may also exhibit different gating changes, such 
CIS db change in the channel activation threshold. 

According to the model presented, the action potential propagates as a result of mechano- 
electrical activation of Nav channels by both straggling axoplasmic ionic currents and an 
axoplasmic pressure wave. Further, the transmembrane voltage spike elevates intracellular 
Ca 2+ concentration ([Ca 2+ ];), either by activating membrane Cav channels and allowing 
entry of external Ca 2+ into the cell, or by Ca 2+ release from intracellular stores. Then, free 
intracellular Ca 2+ initiates a contraction in the filament network anchored to the membrane 
interior. Ca 2+ is a good candidate for transducing the H-H voltage spike into filament 
contraction for several reasons. First, Ca 2+ ion channels are ubiquitous. To quote from 
Ref. [27J, "voltage-gated Ca channels are found in every excitable cell. In fact, I feel they 
define excitable cells." Ca 2+ channels are often overlooked, as they are usually found in 
low density, with typical ionic fluxes during an action potential 100-fold lower than Na 
fluxes {27)]. Furthermore, many types of Ca 2+ ion channels are directly and rapidly gated 
by voltage {28 1. Also, an increase in [Ca 2+ ]; is often associated with initiation of motion in 
cells, from motility in freely moving cells and muscle contraction to synaptic vesicle release at 
synapses. [Ca 2+ ]; in a resting cell is kept extremely low (30-200 11M) by the combined actions 
of a Ca 2+ pump and a Na + -Ca 2+ exchange system on the surface membrane and by Ca 2+ 
pumps on intracellular organelles; extracellular Ca 2+ concentration ([Ca 2+ ] ) is normally 
much higher. Ca 2+ ions that influx during the action potential should instantly increase the 
Ca 2+ level many fold in the vicinity of the filament network attached to the inner membrane 
surface; that, coupled with the ability of Ca 2+ to induce rapid conformational changes in 
proteins (such as upon its binding to actin filaments in the actin-myosin muscle complex), 
could provide a fast contractile response necessary to amplify a propagating pressure pulse. 
Experiments show a drastic change in nerve impulse conduction in low [Ca 2+ ] D , with the 
conduction poor or blocked at initiation of cathodal current; these nerve fibers do respond 
better to anodal break excitation, which provides a higher gradient for Ca 2+ entry 29]. 
However, both [Ca 2+ ]; and [Ca 2+ ] Q significantly affect gating of voltage-gated channels (e.g., 
[Ca 2+ ] meaningfully shifts gating thresholds of all voltage-gated channels) [27|, which affects 
propagation in the pure H-H framework as well. This may be a reason why measurements 
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of nerve impulse conduction velocity with depleted [Ca 2+ ] D are not readily available in the 
literature. We note that, since the radial bilayer contraction follows the pressure pulse at 
the same velocity, even a small contraction that yields a small amplification will accumulate 
over distance and, thus, create a larger amplifying effect. Within this context, it should 
also be noted that a pressure pulse that propagates in an inviscid incompressible fluid that 
is enclosed in a viscoelastic tube will dissipate energy and subsequently decay. This decay 
occurs because the pressure exerted by the tube on the fluid during tube radial expansion is 
greater than that exerted during tube radial contraction. As a result, the total work done by 
the tube on the fluid is negative. In axons, the situation may be the reverse, in which a larger 
pressure exerted by the bilayer during its contraction may result in the overall transmission 
of kinetic energy to the axoplasm. We also posit that generation and amplification of 
the axoplasmic pressure pulse may proceed through electro-mechanical coupling such as 
voltage- induced membrane movement [30] resulting from the H-H voltage spike. Apart from 
membrane movement, it is plausible that elevated [Ca 2+ ]; associated with a propagating 
action potential may act on interior actin-myosin and other cytoskeletal elements to cause 
local cytoplasmic pressure variations that, accumulated over distance, may lead to generation 
of the axoplasmic pressure pulse. 

III. PROPAGATION OF AXOPLASMIC PRESSURE PULSES 

We consider the propagation of small amplitude, axially symmetric pressure pulses in 
a viscous compressible axoplasm enclosed in a circular cylindrical thin-walled distensible 
membrane. The theoretical analysis of waves in viscous fluids enclosed in tubes has been 



presented in Refs. 



31 



U-I33I] . For the limiting case of a viscous compressible axoplasm in a rigid 
tube, the pressure pulse is an axoplasmic density disturbance that propagates along the axon 
similarly to a sound wave packet, with the pulse potential energy stored in the axoplasm 



bulk deformation 



31 



33l | . In the opposing limit of a viscous incompressible axoplasm in a 
distensible membrane, the pressure pulse manifests itself as an increase in the axon diameter 
and associated membrane area expansion, with the pulse potential energy stored in the 
membrane strain 32, 3^]. Let us consider a pressure pulse of central frequency u propagating 
in the axoplasm of density p, compressibility k and dynamic viscosity fi, enclosed in a 
cylindrical membrane of radius R, thickness h and area expansion modulus K. Let us also 



7 



assume that for the deformation of area expansion (when the membrane material behaves 
in an elastic solid manner), the membrane Young's modulus is E and the Poisson's ratio is 
v. If viscosity is neglected, the pulse propagates with the velocity [34j |: 



v 



(i) 



For typical unmyelinated axons, i.e., assuming K = 0.8 N/m (twice the single plasma 
membrane area expansion modulus of 0.4 N/m 35], to account for the combination of the 
axon and adjoining glial cell membranes), k = 4.04 ■ 10~ 10 Pa -1 (water at 38 °C l^). 
p = 1050 kg/m 3 (axoplasm jsjzj]), and 2R = 1 //m, Eq. (1) yields the pressure pulse velocity 
of 28 m/s. Note that if the axon membrane is assumed to be indistensible (^P <C k), Eq. 
(1) yields a pulse velocity of 1535 m/s. 



A. Viscous axoplasm 



Viscous forces in the axoplasm and in the axon membrane in general decrease the velocity 



of axoplasmic pressure waves and introduce a decay [31H33J. The membrane viscous tension 
is characterized by the ratio of the surface viscosity (on the order of 10 -6 N-s/m 38]) and the 
time scale of the pressure pulse (here assumed to be close to the action potential duration, 
about 1 ms), yielding 10~ 3 N/m. Since this value is much smaller than the membrane 
elastic area expansion modulus, 0.4 N/m, the effects of membrane viscosity on the pulse 
propagation can be neglected 1 



Concerning the viscosity of the axoplasm, for the general 
case of a viscous compressible fluid enclosed in an elastic tube, the dispersion equation is 
quite complicated and has not been solved in closed form. However, the equation is simplified 



and can be solved in the high viscosity limit, defined by 3lM33]: 



a: 



rJ^ « 1. 
V ^ 



(2) 



Here we evaluate Eq. (2) for typical axons. As before, we assume that the duration of the 
axoplasmic pressure pulse is similar to that of the action potential, about 1 ms, yielding 
u ~ 3142 rad/s. Axoplasm viscosity p, as is relevant for perturbations associated with 
small-amplitude axoplasmic pressure pulses, is assumed to be similar to the viscosity of 
water, p = 6.82 • 10~ 4 Pa-s at 38 °C 36]. Taking 2R = 10pm (the typical internal myelin 
diameter for Aa myelinated fibers) and p = 1050 kg/m 3 yields a = 0.35, which should be 
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small enough for high-viscosity approximation to the dispersion equation to hold [33J for 
axons of similar or smaller diameters. When a ^> 1 (2R ^> 29/im), the effects of viscosity 
are small, and the pressure pulse propagates with little decay and with the velocity described 
by Eq. (1). 

When Eq. (2) holds, the phase velocity of axoplasmic pressure waves (the velocity of the 



harmonic waves) is given by 3lM33|: 



Vph 



•at-Vo, 



(3) 



where a and vq are given by Eqs. (2) and (1), respectively; c 



1 when ^ < k (the 



limit of a rigid membrane), while c 



when ^ 3> k (the limit of a soft membrane). 



v / 5 = 4^ """" K 

The factor c accounts for movement of the membrane in the axial direction caused by the 
viscous axoplasm in the limit of a soft membrane; here and below we assume v = 0.5 for the 



incompressible lipid bilayer, leading to 



1.15. The velocity of the pressure pulse (the 



"group" velocity) is given by v gr = |£, where k is the wavenumber. From Eqs. (1-3), and 



using the relations v p h = f and v gr 



dm 

dk ■ 



we obtain the velocity of the axoplasmic pressure 



pulse in the high viscosity limit: 



C ■ a ■ Vq = c ■ R 



(4) 



Under the same 
decay length), is 



imit, the distance over which the pulse amplitude decreases e-fold (the 
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-|33j: 



(5) 



B. Myelinated axons 



In the H-H model of propagation of the action potential, the role of the myelin sheath is to 
insulate the membrane electrically and decrease its electrical capacitance; both these effects 
contribute to more rapid conduction of the action potential. However, the myelin sheath 
is also quite rigid. About 80% of dry myelin content is lipid, with a roughly 2:2:1 molar 



ratio of the three major lipids, cholesterol, phospholipids and galactolipids [39]. Cholesterol 
is known to increase the strength and elastic modulus of lipid bilayer s, with the highest 
modulus and highest strength achieved for 50 mol% cholesterol 40], i.e., rather close to 



the concentration in myelin. The remaining 20% of dry myelin content consists of protein, 



including various types of fibrous proteins that can form rigid structures. Altogether, for the 
deformation of circumferential strain, the myelin elastic area expansion modulus is about 2 
N/m for a single myelin layer of 4 nm thickness j^oj] . Noting that K = Eh (for v = 0.5), the 
myelin effective Young's modulus is E = 5 ■ 10 8 Pa. The ratio of the axon diameter (internal 
in diameter) to the nerve fiber diameter (external myelin diameter) is typically about 



mye 
0.7 



4l|, implying a myelin thickness of h f» 0.43.R. Therefore, =^ = W- ~ 9.3 -10 9 Pa 1 . 
This is about 23 times greater than the compressibility of water at 38°C, k = 4.04 • 10~ 10 



Pa -1 . Hence, t>o, as used in Eq. (4), is only a/24 ps 4.9 times less than the speed of sound 
in water (see Eq. (1)). In other words, the myelin sheath is so rigid that its distensibility 
decreases the pulse speed only about 4.9 times compared to the theoretical maximum if axon 
walls were absolutely rigid. For typical myelinated axons, e.g., cat myelinated axons of 10 
/iin fiber diameter at 38 °C, Eq. (4) yields v gr = 88 m/s (using E = 5 ■ 10 8 Pa, 2R = 7/iin, 
h = 0A3R, v = 0.5, co = 3142 rad/s, // = 6.82 • 10~ 4 Pa-s and k = 4.04 • 10~ 10 Pa" 1 ). Given 
the roughness of the estimate, this theoretical value of 88 m/s is in excellent agreement with 
the measured action potential velocity of 60 m/s for these fibers 42j. From Eq. (5), the 
decay length for the propagation of the pulses in myelinated segments of such a fiber (where 
the pulse propagates passively) is 14 mm. Assuming that the distance between the nodes 
of Ranvier is 1 mm (100 times the external diameter of myelin |43]), the pulse amplitude 
decays 7% over the internodal distance. 



C. Unmyelinated axons 



Unmyelinated nerve fibers usually have diameters of 0.1 - 1.2 /im 44J. The action po- 
tential conduction velocity, in m/s, is given approximately by v ~ 1800a/R, where R is the 
axon radius in meters 44| (values quoted in the literature range from v ~ 1000-\/i? 45| to 



3000V.R 



46||). Therefore, for an axon of 1 /im diameter, conduction velocity is about 



1.27 m/s. Assuming 2R = 1 /im, K = 0.8 N/m (twice the single plasma membrane value of 
0.4 N/m 35|), and u, [i, k and v as before for myelinated axons, Eq. (4) yields a pressure 
pulse velocity of v gr = 1.11 m/s, very close to the experimental value. From Eq. (5), the 
decay length for this fiber is 0.18 mm. It is truly remarkable that, for both myelinated and 
unmyelinated fibers, the predicted velocity of propagation of axoplasmic pressure pulses is 
very close to the measured velocities of nerve impulses. 
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D. Dependence of propagation velocity on fiber diameter 



As can be seen from the above numerical estimates, for both myelinated and unmyelinated 
axons, the condition ^ 3> k is true; and therefore, Eq. (4) simplifies to: 



This relationship shows that for unmyelinated axons, the axoplasmic pressure pulse velocity 
is proportional to the half-power of the axon diameter, assuming that the pulse duration, 
the membrane area expansion modulus and the axoplasm viscosity do not vary with axon 
diameter. For myelinated axons, introducing 7 = h/R (the ratio of the myelin sheath 
thickness to the axon radius) and using K = Eh, we can rewrite Eq. (6) as: 



v 



! 2cuE~f 



Therefore, for myelinated axons, assuming that pulse duration and specific properties of 
axoplasm and myelin are independent of diameter and that the myelin sheath thickness 
scales in proportion to the fiber diameter (i.e., 7 is constant), the pressure pulse velocity 
scales linearly with the fiber diameter. It is interesting that a linear dependence of pressure 
pulse velocity on fiber diameter is also obtained for myelinated axons without assuming 
K < ^, from Eq. (4): 



V = c ■ R U 2 . (8) 

In summary, given that specific properties of axoplasm and myelin are independent of fiber 
diameter and that the myelin sheath thickness scales in proportion to fiber diameter, the 
axoplasmic pressure pulse velocity scales as yR for unmyelinated axons and as R for myeli- 



nated axons, similar to what is measured experimentally 



4J-|46| 



E. Dependence of propagation velocity on temperature 

The duration of the action potential increases with decreasing temperature, e.g., with Q 10 
(the ratio of the duration at one temperature to the value at a temperature 10 °C lower) 
of 3.4 at 37 °C for cat vagus myelinated fibers [471 ] . Assuming that the duration of the 
axoplasmic pressure pulse is similar to that of the action potential, taking the change in /1 
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from 37 °C to 27 °C as for water, Qio = 0.81 [36], and assuming that the other parameters 
in Eq. (4) do not change with temperature, we obtain the predicted Qio of 2 for the nerve 
impulse propagation velocity. This is rather close to the experimental value of 1.6 



48]. 



IV. MEYER-OVERTON RULE AND THE ACTION OF ANESTHETICS 



If axoplasmic pressure waves are involved in a meaningful way in a normal propagation of 



the nerve impu 



anesthetics 



16 



se, we arrive at a possible explanation of the Meyer-Overton rule for inhaled 



a-ila]. Essentially, the rule states that anesthetic potency is determined b y th e 
bilayer concentration of anesthetic, independent of its molecular identity (e.g., Ref. [24j). 
Therefore, the Meyer-Overton rule can be explained if there is a property of the lipid bilayer 
that is essential for nerve impulse propagation and that, in addition, depends only on the 
bilayer concentration of anesthetic. As is shown below, one such property may be the 
bilayer area expansion modulus, which affects both the decay and the velocity of axoplasmic 
pressure pulses. It is also shown that anesthetics may inhibit the mechanical modulation 
of Nav channels by the traveling pressure pulse (and therefore contribute to suppression of 
conduction of the action potential) in a manner consistent with the Meyer-Overton rule. 

Suppose that an anesthetic molecule is dissolved within the lipid bilayer (anywhere from 
the headgroup to the hydrophobic region), and it weakens and destabilizes the bilayer struc- 
ture nonspecifically, e.g., through creating a defect in the packing of lipid molecules and 
thereby nonspecifically disordering lipid chains (Fig. 4(a)). Such disordering interaction 
will generally decrease the bilayer elastic area expansion modulus in the vicinity of the 



lipid/anesthetic molecular complex (it is well established that bilayers with 
chain disordering generally display lower elastic moduli and tensile strengths 



arger lipid 



40j). Using a 



simple macroscopic model of the membrane 40j, shown in Fig. 4(b), the total membrane 
area expansion modulus can be expressed as: 

fi+i^y 1 , (9) 

-"■M K L/A/ 

where the right-hand side of the equation is a combination of the area expansion moduli 
of the two components, namely, non-anesthetized membrane Km (which describes the elas- 
ticity of the lipid bilayer, embedded proteins and other native structures) and a purported 
lipid/anesthetic molecular complex Kl/a (Kl/a < Km), scaled by their area fractions clm 
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FIG. 4: (a) Disordering of lipid packing order by an anesthetic molecule dissolved in the bilayer. 
(b) Pockets of altered area expansion modulus in the vicinity of anesthetic molecules. 

and cll/a- Assuming that anesthetic molecules are dissolved only in the lipid segments of 
the membrane and that they do not cluster together, the area fractions au and a L / A are 
given by: 

ah/A = 1 - a M = n A ■ h ■ A L/A , (10) 

where n A represents the concentration of anesthetic molecules within the membrane, in units 
of molecules per volume, h the bilayer thickness, and A L / A the membrane area occupied by a 
single lipid/anesthetic molecular complex. Note that Al/ A and Ki/ A are independent of the 
chemical identity of the anesthetic, due to the purported nonspecificity of anesthetic/lipid 
interaction. Combining Eqs. (9) and (10): 

„ (I - n A hA L/A n A hA L/A \~^ 

K = { Km + -l^T) • (U) 
13 



where the right-hand side of the equation does not depend on the identity of the anesthetic, 
but only on its bilayer concentration n\- As was shown above, both the velocity and the 
decay length of axoplasmic pressure pulses vary as \J~K. Therefore, a decrease in the area 
expansion modulus K due to an increase in the anesthetic concentration ha will reduce the 
pressure pulse velocity and increase its decay, which, when ua is of sufficient magnitude, may 
lead to an inhibition of the pulse propagation. From the form of Eq. (11), it follows that 
anesthetics acting in such a way will comply with the Meyer- Overton rule. From the line of 
argument above, it is clear that any nonspecific anesthetic/lipid interaction that results in 
decreasing K with increasing via and thus leads to a block in nerve impulse conduction will 
be consistent with the Meyer-Overton rule. 

Another possible mechanism of action of anesthetics consistent with the Meyer-Overton 
rule relies on inhibition of modulation of Nav channels by the traveling pressure pulse, in 
cases where such modulation is mediated by changes in the membrane lateral pressures 
(Fig. 3(a)). As before, we assume that anesthetic molecules dissolved in the bilayer weaken 
it nonspecifically and that increasing anesthetic concentration ua decreases the membrane 
area expansion modulus K. As was discussed above, the potential energy of a propagating 
axoplasmic pressure pulse is mainly stored in the membrane strain. The membrane strain 
energy, E p , is proportional to LAt 2 /K, where L is the pulse length and At is the peak in- 
crease in the membrane tension r from equilibrium. Since L ~ \[K (assuming a fixed pulse 
duration, see Eq. (6)), for the same pulse energy (E p fixed), it follows that At ~ K 1 ^ 4 . 
Therefore, a decrease in the area expansion modulus K will result in a reduced peak tension 
increase At. At sufficiently high ha (low K), the decreased At may not suffice to stretch- 
modulate Nav channels to the degree that is required for the process of co-propagation of 
the electrical and mechanical excitations to occur, possibly leading to conduction block. It 
should also be noted that, if the anesthetic/lipid interaction changes the lipid bilayer area 
without affecting its area expansion modulus, the equilibrium tension in the membrane may 
change, which, in turn, should affect gating of Nav channels and also possibly lead to conduc- 
tion block. This anesthetic/lipid interaction would also conform to the Meyer-Overton rule 
if it is nonspecific. In this context, it should be noted that an anesthetic/lipid interaction 
that alters a mechanical property of the lipid bilayer in a nonspecific manner would probably 
involve a defect-like structure, in which physical properties of the anesthetic compound are 
not important provided that the compound satisfies certain criteria for creating the defect 
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(e.g., being hydrophobic, nonpolar, and nonreactive and possessing certain spatial and mass 
dimensions) . 



V. DISCUSSION 



It has recently been suggested 



0,Q 



in 



that the H-H model [l[ of the action potential 



does not provide a satisfactory description of the nerve impulse because it does not include 
the mechanical {4-8] and optical Q, Q changes associated with the action potential and 
because it is inconsistent with observed reversible changes in temperature and heat [q 13 ] . 



Furthermore, the famous Meyer-Overton rule 16Ml8l| is incongruous with protein models of 

J li 

action of anesthetics [3j, which are consistent with the H-H model. A theory of the action 
potential based on mechanical solitons propagating in the axon membrane was suggested 
instead 3J]. Although this theory is in principle able to account for the above-mentioned 
phenomena, apparently it cannot explain changes in the propagation velocity that are related 
to axon diameter. 

Following related earlier work 



49 



50j, here we propose that the mechanical, optical 
and thermal changes associated with the action potential 4Ml5| are the result of an ax- 
oplasmic pressure pulse that accompanies the action potential and has a roughly similar 
duration. We have shown that, for both myelinated and unmyelinated axons of different 
diameters, such pressure pulses propagate with a velocity close to the measured action po- 
tential velocity. Hence, any axoplasmic disturbance resulting from a propagating action 
potential (e.g., caused by the influx of extracellular Ca 2+ ions and their presumed action on 
acto-myosi n cy toskeletal elements, or the disturbance caused by volt age- induced membrane 



movement 



301]). should accumulate into a larger, shock-like axoplasmic pressure wave such 



as we propose. Further, based on known axon properties, we suggest a model in which the 
axoplasmic pressure wave accelerates propagation of the H-H type excitation wave, leading 
to co-propagation of the two waves at the velocity of the axoplasmic pressure wave. We 
also suggest several mechanisms for pressure wave amplification that could counteract its 
expected viscous decay. The presented model has experimental basis in the above-mentioned 
observations of mechanical, optical and thermal changes associated with the action poten- 
tial, which all point to the involvement of a mechanical or thermodynamic process in the 
propagation of the nerve impulse. Importantly, the simple and powerful Meyer-Overton rule 



15 



16 



18] , valid over five orders of magnitude, currently appears to be difficult to interpret in 



the context of presently accepted models. Within the presented framework, however, the 
Meyer-Overton rule can be explained based on known interaction of anesthetics with lipid 
bilayers. A simple test of the presented theory is its ability to predict the velocity of the 
nerve impulse. We have shown that, over two orders of magnitude (for myelinated and un- 
myelinated axons), the predicted velocity is close to that which is experimentally observed, 
despite the fact that the precision of calculations is limited by current knowledge of the 
mechanical properties of axons. The dependence of nerve impulse velocity on fiber diameter 
and temperature is also reproduced well by the model. 

If some types of neurons or other excitable cells indeed rely to a significant degree on 
the presented mechanism of signal transmission, perhaps so much so that propagation of 
the signal within the pure H-H framework is not feasible, it becomes plausible that various 
pressure-mediated processes may play important roles in the function of these cells. De- 
polarization of the axon hillock to the excitation level, which presumably leads to a radial 
contraction of the hillock membrane and the initiation of an axoplasmic pressure pulse, may 
proceed not only through a passive spread of ionic currents from the postsynaptic sites, but 
also through mechanical modulation of gating of ion channels in the vicinity of the axon 
hillock. Therefore, synaptic integration may include relay of cytoplasmic pressurization and 
displacement from postsynaptic sites through the dendritic tree to the axon hillock. It follows 
that learning and memory in neurons may be partially reflected in the strength and duration 
of a pressure-generating contraction of a protein filament network that is anchored to the 
postsynaptic membrane; the contraction may be activated as a result of the development 
of postsynaptic potentials. Hence, some of the structural elements that underlie learning 
and memory may be similar to those that regulate muscle cell strength (e.g., the number of 
actin-myosin filaments), which would provide a possible link between learning and memory 
in neurons and the changes that occur in muscle cell strength due to exercise. The existence 
of a dense undercoat ing just beneath the axon membrane in the axon hillock and nodes 



of Ranvier (e.g., Ref. (5l|), currently functionally unexplained, may in fact be the putative 
contractile filament network attached to the inner membrane surface, expressed specifically 
in regions where energetic generation and amplification of the axoplasmic pressure pulses 
has to occur. Nerve impulses generated in the axon hillock may affect processes in the 
neuronal cell body, affecting dendrites and contacting synapses through back-propagating 
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cytoplasmic pressure pulses. Action potential generation and propagation may be affected 
as a result of stretch-modulation of membrane voltage-gated channels by bilayer stresses 
resulting from altered [Ca 2+ ]; levels and their action on cytoskeletal elements following pre- 
vious spike activity. Pressure pulses associated with action potentials may play a role in 
transport of substances in excitable cells. Synaptic vesicle release at axon terminals may be 
facilitated by a contraction of the knob-shaped membrane skeletal protein network of the 
synapse, creating a force pushing the vesicles toward the synaptic cleft. The major compo- 
nents of the model, propagating hydraulic pulses and radially contracting tubes generating 
an axial force on the tube contents, are exploited in the cardiovascular and gastrointestinal 
systems, suggesting a possible evolutionary link between these systems and excitable cells. 
Finally, it is suggested that the evolutionary predecessor of the ability to transmit nerve im- 
pulses may be the well-known ability of freely moving cells to control membrane movement 
by forces of microfilaments located just beneath the plasma membrane (e.g., Ref. {oil), in 
response to external as well as internal mechanical stimuli. 
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